arXiv: 1507.05514vl [hep-ph] 20Jul2015 


JUNO and Neutrinoless Double Beta Decay 

Shao-Feng Ge* 1 and Werner R.odejohanrd 1 


^Max-Planck-Institut fur Kernphysik, Postfach 103980, 69029 Heidelberg, Germany 

July 21, 2015 


Abstract 

We study the impact of the precision determination of oscillation parameters in the JUNO experiment on half-life 
predictions for neutrinoless double beta decay. We show that the solar neutrino mixing angle can be measured by 
JUNO with below 1% uncertainty. This implies in particular that the minimal value of the effective mass in the 
inverted mass ordering will be known essentially without uncertainty. We demonstrate that this reduces the range 
of half-life predictions in order to test this value by a factor of two. The remaining uncertainty is caused by nuclear 
matrix elements. This has important consequences for future double beta decay experiments that aim at ruling out 
the inverted mass ordering or the Majorana nature of neutrinos. 


1. Introduction 


With neutrino mass and lepton mixing firmly established as facts, a minimal basis to describe those phenomena has 
been developed. In this “3 Majorana neutrino paradigm”, three massive Majorana neutrinos are described by the 
neutrino mass matrix m„: 



( 1 . 1 ) 


There are nine physical parameters in to„, usually parametrized as 3 masses, 3 mixing angles and 3 phases [1]. Within 
this simple framework, 3 tasks are eminent: determining the parameters as precisely as possible; checking if the 
minimal description (3 Majorana neutrino paradigm) is correct; explaining the measured parameter values. 

Yet to be determined experimentally are the neutrino mass ordering, the octant of the atmospheric mixing angle, 
the Dirac CP phase, the absolute mass scale, whether neutrinos are Dirac or Majorana type fermions, and the 
Majorana CP phases if neutrinos are Majorana type. The first three play a role in neutrino oscillations, with which 
one can determine in total 6 parameters, 2 mass-squared differences, 3 angles and the Dirac phase. For the remaining 
parameters and properties non-oscillation experiments are inevitable. While our current information on the established 
oscillation parameters is already impressive [2], further improvement is of course needed. For instance, this is necessary 
to rule out flavor symmetry models (see e.g. [3]) or to check for new physics such as non-standard interactions, unitarity 
violation, long-range forces, etc. In this paper we address however another aspect of precision determination of neutrino 
oscillation parameters, namely its impact on neutrinoless double beta decay (0^2/3) [4, 5]. 

It is well-known that 0v2/3 can in principle contribute to determining the Majorana or Dirac nature of neutrinos 
and also to the question of the neutrino mass ordering: in the standard paradigm the effective mass, on which the 
amplitude of the process depends, cannot vanish if the inverted mass ordering is realized. Therefore, if an oscillation 
experiment shows that neutrino masses have an inverted ordering, we know that the process must happen with a 
certain half-life Not observing the decay with a half-life limit above means that neutrinos are Dirac 

particles. In turn, not knowing the mass ordering and not observing the decay with a half-life limit above rules 

out the inverted ordering in case neutrinos are Majorana particles 1 . The timescales of determining the mass ordering 
[6] and reaching half-life limits in the inverted ordering regime are comparable [ ] and subject to uncertainties, so 
both possibilities could happen. The effective mass value that future experiments need to reach in order to fulfill the 
two goals mentioned above needs to be known as precisely as possible, since it enters the half-life quadratically. If the 
experiment is dominated by background the situation is even worse. 
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x Of course, new physics in neutrinoless double beta decay such as light sterile neutrinos or TeV-scale left-right symmetric models could 
modify such statements, see Refs. [ , ] for a general discussion. We assume here that no contribution other than three massive Majorana 
neutrinos to double beta decay plays a role, thus we stay within the best motivated (standard) interpretation of the decay. 
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As argued in [8], the upper limit on the minimal half-life in the inverted ordering depends strongly on the solar 
neutrino mixing angle and its current uncertainty introduces a sizable uncertainty in half-life predictions. This “particle 
physics uncertainty” is of the same order as the “nuclear physics uncertainty”, i.e. the notorious nuclear matrix 
elements. Both uncertainties should be reduced. In this paper we demonstrate that future precision data on the 
oscillation parameters essentially removes the particle physics uncertainty. Towards this end, we apply the NuPro 
] package, written by one of us (SFG), to evaluate the precision with which the upcoming JUNO experiment [10] 
will determine the neutrino oscillation parameters and thus in particular the minimal effective mass in the inverted 
ordering. For demonstration, we will use JUNO as an example, we note that RENO-50 [11], with a very similar 
configuration, will reach similar precision. 

We demonstrate that such medium baseline reactor experiments can measure the solar neutrino mixing angle with 
unprecedented precision of even less than 1 %, in addition to their main goal of measuring the neutrino mass ordering. 
We translate this precision in half-life ranges of neutrinoless double beta decay for various isotopes and nuclear matrix 
element calculations. This will help 0v2(3 experiments to possibly rule out the Majorana nature and evaluate their 
requirements to achieve this goal. 

JUNO can not only determine the mass ordering, but also can measure the solar neutrino mixing angle with 
impressive precision at the same time. Thus, a single experiment can provide double beta decay experiments with 
all necessary information for ruling out the Majorana nature. To be more precise, the minimal value of the effective 
mass in the inverted mass ordering depends on the atmospheric mass-squared difference (very weakly on the solar one 
as well), 0 i 2 and $ 13 , i.e. the same set of parameters the electron neutrino survival probability in reactor experiments 
depends on. Thus, the latter have a direct correspondence to neutrinoless double beta decay experiments (and in 
principle also to single beta decay experiments). 

Of course, if JUNO determines that neutrinos enjoy a normal mass ordering, the interesting link to double beta 
decay is somewhat lost, though the uncertainty of the effective mass in the normal ordering will still be reduced 
significantly. In case JUNO would not be able to determine the mass ordering because of limited energy resolution, 
we show that the precision on the solar neutrino mixing angle is not affected. 

This paper is organized as follows. In Sec. 2 we discuss 0u2P with three Majorana neutrinos and how the ability of 
ruling out the inverted ordering or the Majorana nature is related to the uncertainty of neutrino oscillation parameters. 
Then, we evaluate the precision that can be achieved at JUNO in Sec. 3, which is used in Sec. 4 to derive the required 
half-life sensitivity that C02/3 experiments need to provide. Our conclusion can be found in Sec. 5. 


2. Effective Neutrino Mass and Half-Life 

Neutrinoless Double Beta Decay (C02/3) is fundamentally important to particle physics as its observation implies 
lepton number violation, similar to the baryon number violation implied by an observation of proton decay. Several 
mechanisms for 0u2 /3 exist [4, 5]. In the best motivated interpretation, C02/3 is mediated by Majorana neutrinos with 
mixing observed in neutrino oscillation experiments. The decay half-life can be expressed as 

(Ti 0 ; 2 ) _1 = G 0v \M 0v \ 2 (m v } 2 , (2.1) 

where G°" is the well-known phase space factor and M 0v the nuclear matrix element. Neutrino mass and mixing 
enters through the effective electron neutrino mass ( m v ) 

(m v ) = \c 2 s c 2 r mi + s 2 s c 2 r m 2 e la + s 2 mze' ,p \ , (2.2) 

with ( c a ,s a ) = (cos0 Q ,sin0 a ). Here, we have adopted the notation for mixing angles according to 9 a = 023 denoting 
the atmospheric mixing angle, 6 r = 0 i 3 the reactor mixing angle and 0 S = 0 i 2 the solar mixing angle. For (02/3, the only 
relevant mixing angles are 0 S and 9 r . Apart from these, the relevant parameters from the particle physics side involve 
the three neutrino mass eigenvalues rrii and the two Majorana CP phases a and /3. Measured by neutrino oscillation 
experiments are the mixing angles and the two mass-squared differences, Am 2 = m 2 — m 2 and Am 2 = I 7713 — m 2 \. 
Note that for the mass-squared difference between m\ and to| only the magnitude has been measured. The different 
sign of Am 2 leads to quite different mass patterns, the normal ordering (NO) with m\ < m 2 < to 3 and the inverted 
ordering (10) with m 3 < m\ < m 2 . Fig. 1 shows for both mass orderings the effective mass versus the smallest mass, 
as well as versus the neutrino mass parameters that are accessible in direct searches and cosmology. Both the current 
as well as future 3<r ranges, to be determined in the later part of the paper, are given. 

Even with the absolute mass scale and mass ordering fixed, the effective electron neutrino mass (m v ) can still vary 
as a function of the two unknown Majorana CP phases a and /3. For IO, the maximal and minimal values are 

(mv )max = \Jml + A m 2 a c 2 s c 2 r + yjmj + Am 2 a + Am 2 s s 2 c 2 + m 3 s 2 , (2.3a) 
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Fig. 1: The effective mass as a function of the smallest mass eigenvalue (mi for NO and m 3 for 10), /3-decay mass 
nif 3 = VWi 1 2 m|, and the sum of mass eigenvalues, within the 3er range before (prior) and after (posterior) JUNO. 

(w„)n° n = \J m| + Am 2 c 2 c 2 - yjml + Am 2 + Am 2 s 2 c 2 - m 3 s 2 . (2.3b) 

Since m 3 is the smallest mass, s 2 <C c 2 , c 2 ~ 2s 2 and A m 2 Am 2 , the first term actually dominates over the second. 
As the two limits in (2.3) increase with m 3 , there is an universal lower limit [ 12 ] for 10 when m 3 approaches zero (the 
lower limit (m„)^ n can also easily be extracted within a geometrical picture [13]): 

("C>Sn V Am l( c2 s ~ s 2 s) C r ■ (2-4) 

This value can put an upper bound on the half-life time T®j 2 through (2.1). To exclude the inverted ordering, it 
is necessary for C02/3 experiments to go beyond this limit. How precisely do we know this limit? Currently, the 
atmospheric mass-squared difference A m 2 and the reactor mixing angle O r have been measured with precision at the 
few percent level. However, the solar mixing angle 9 S can contribute a significant uncertainty to (m y ) [8]. 

According to the latest global fit [2], the 3<r range of the solar mixing angle is 0.444 < cos 2 6 S < 0.250. This implies 
at 3cr C.L. a total uncertainty of 

= (0.0127... 0.0198) eV, (2.5) 
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Fig. 2: Current limits on the effective mass (90% C.L.) for some isotopes for the range of nuclear matrix elements 
given in Table 1. The vertical lines are current experimental half-life limits. 
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Fig. 3: Experimental constraints (90% C.L.) on the effective mass for some isotopes for the range of nuclear matrix 
elements given in Table 1. The vertical lines are current experimental half-life limits, the horizontal band is the current 
range for the minimal value in the inverted ordering, the lower solid line is after JUNO. The upper solid line is for the 
half-life to enter the inverted ordering regime, = sj Am^. 


which varies by a factor of 1 . 6 . Apart from (to,,), there is of course another source of uncertainty in predicting Tj 0 /, 
namely the nuclear matrix elements . Their uncertainty, see Table 1, is of similar magnitude as the one coming 
from the oscillation parameters, in particular 9 S . Nuclear physics uncertainties will be addressed by the nuclear physics 
community. Here we focus on the particle physics uncertainty in (m„), where progress is essentially guaranteed. In 
Fig. 2 we show, using the range of nuclear matrix elements from Table 1, for a subset of isotopes the current effective 
mass limits. The limits are 

(to„) < (0.22...0.53) eV 
(mj) < (0.36...0.90) eV 
(mj) < (0.27... 1.00) eV 
(mj) < (0.15...0.35) eV 

see also [ ]. Note that these limits are not just for the three light neutrinos from Eq. (1.1), but include in principle 
possible sterile neutrinos up to masses slightly below the Fermi scale, i.e. the constraint is actually on )T) 

1,71 

where the sum goes over all neutrinos below m n ~ 50 MeV. Note finally that the matrix elements depend roughly 
quadratically on the axial coupling constant gA , and thus the effective mass limits in Eq. (2.6) would weaken by 
(1.27/g a) 2 in case gA would quench, i.e. become smaller as the nuclear mass becomes larger [15], see also [16]. The 
half-life limits we obtain in what follows would become stronger by a factor (1.27/< 7 ^ 4 ) 4 . Fig. 3 shows for a larger 
subset of isotopes the effective mass limits as a function of half-life limits, showing in particular the necessary half-lifes 
to rule out the inverted ordering, to be derived in Sec. 4. We also give the necessary half-lifes to enter the regime in 
which for the inverted ordering in case of vanishing smallest mass the effective mass should lie, (TO„)™ ax = \J Ato 2 c)(. 

Note that the effective mass enters the half-life quadratically. If a 0v2/3 experiment is dominated by background, 
the half-life that an experiment can reach is proportional to [17] 


for 76 Ge, 
for 100 Mo, 
for 130 Te, 
for 136 Xe, 


( 2 . 6 ) 


Xj°/ 2 oc a x e x 


M xt 
B x AE 


(2.7) 


where a is the isotopic abundance of the double beta emitter, M the fiducial mass, t the measuring time, e the 
detection efficiency, A E the energy resolution around the peak, and B the constant background rate in units of 
counts/keV/kg/yr. We see that a factor of 2 uncertainty in the effective mass implies experimentally either a factor 
4 in half-life or a combined factor of 16 in the experimental parameters in the square root in (2.7). For the current 
3cr range, a variation by a factor of 1.6 in (TO„)[9 n corresponds to a factor of 2.56 variation in the half-life, and thus a 
factor of 6 in the experimental parameters. Obviously this is challenging, especially when it comes to estimating the 
necessary size of the detector and the runtime to reach the limit of excluding IO. 
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In the next Section we will demonstrate how the precision on the effective mass can be improved at JUNO. 


3. Precision Measurement of Solar Mixing Angle at JUNO 


Reactor neutrinos are generated by nuclear reactions. Hence, their typical energy is in the MeV range which is not 
enough to produce muon or tau leptons in the final state. The only accessible oscillation channel is the electron 
survival probability P ee , which depends on only two mixing angles, 9 r and 9 S , as well as two mass-squared differences: 

P ee = 1 — 4c/c 2 s s 2 sin 2 A 21 — 4 c 2 c 2 s 2 sin 2 A 31 - 4s 2 c 2 s 2 sin 2 A 32 , (3.1) 

where A^ = (m 2 — m 2 )L/4E v . Note that, just as for the effective mass, the atmospheric mixing angle 9 a and the Dirac 
CP phase Sd are not involved. At short-baseline reactor experiments the dominant oscillation comes from the last 
two terms which are modulated by the large mass-squared difference Am | 2 ~ A?n 31 = Am 2 . The amplitude ~ c 2 s 2 
is essentially independent of the solar mixing angle 9 S , hence short-baseline experiments can be used to measure the 
reactor mixing angle 9 r precisely. This has been and will be done by Daya-Bay, Double Chooz and RENO. 

On the other hand, at longer baseline the contribution from the second term in (3.1) becomes dominant. On top 
of this slow oscillation, many fast oscillations occur. Note that the last two terms in (3.1) are not exactly the same. 
The difference between A 31 and A 32 is roughly Am 2 /Am 2 ~ 3%. In other words, there are two fast oscillation modes 
with slight difference. This can be used to determine the neutrino mass ordering [18, 19]. Depending on the neutrino 
mass ordering one frequency is faster than the other. To make the picture clear, we can formulate the oscillation (3.1) 
in terms of Am 2 and Am 2 , 

Pee = 1 ^ 4 c/c 2 s 2 sin 2 A 2 1 - 4c 2 s 2 sin 2 |A 3 i | 

- 4s 2 c 2 s 2 sin 2 A 2 i cos ( 2 |A 3i |) ± 2s 2 c 2 s 2 sin( 2 A 2 i) sin ( 2 |A 3i |), (3.2) 


where ± corresponds to NO (+) and 10 (—) respectively [20]. The difference in oscillation probability is equivalent to a 
relative phase shift between the two fast frequencies. Since the relative phase difference is only about Am 2 /Am 2 ~ 3%, 
it is necessary for the energy resolution to be better than 3%. By measuring the relative phase shift between the two fast 
frequencies, the neutrino mass ordering can be determined [21, 22, 23]. This corresponds to measuring the atmospheric 
mass-squared difference Am 2 with precision around Am 2 /Am 2 ~ 3%, so that the tiny difference between Am 31 and 
Am 32 can be seen. To illustrate the picture, we show the event rate in Fig. 4. In the current study we focus on the 
JUNO configuration [10] with a 20kt liquid scintillator detector 52 km away from two reactor complexes with 36 GW 
total thermal power. The effective runtime is roughly 4.9 years (normal run time of 6 years with 300 effective days 
per year) and a detector energy resolution 3 %/^E (MeV) as benchmark. 

One side product from this medium baseline reactor neutrino experiment is that the low frequency oscillation due 
to the second term in (3.1), which is essentially modulated by 4c 2 s 2 , can be used to measure the solar mixing angle 
9 S precisely [24, 20]. This constraint mainly comes from the total event rate as shown in Fig. 4. Consequently, it 
will not be affected much by energy resolution, in contrary to the determination of the neutrino mass ordering. A 
rough estimate gives 1 x 10 5 events to be observed, with a statistically fluctuation of 3 x 10 -3 . The actual uncertainty 
should be even further suppressed since not only the total event rate is measured, but a whole spectrum across the 
whole energy range. The position of the main peak in Fig. 4 can be used to constrain the corresponding oscillation 
frequency Am 2 . Neither of them is sensitive to the energy resolution that the detector can finally achieve. Therefore 
a precise measurement of the solar mixing angle 9 S and the solar oscillation frequency Am 2 is guaranteed. In Fig. 4, 
we also show the dependence on the solar mixing angle by plotting the event rates with central value sin 2 9 S = 0.323 
and la upper limit sin 2 9 S = 0.339, respectively. 

We note that there is a study showing that the 5 MeV reactor anomaly can introduce a small shift in the best-fit 
of 9 S [25]. This anomaly is a bump in the ratio of measured and observed neutrino fluxes around energies of 5 MeV, 
as observed at Daya Bay [26], Double CHOOZ [27], and RENO [28, 29]. For a measurement of 9 S the anomaly can 
change the total event rate which is mainly controlled by 4c 2 s 2 as we explained above. Consequently, the best fit of 9 S 
can be shifted away from its true value, although its uncertainty would not be affected much. This difficulty could in 
principle be overcome by an extra near detector. By comparing the event rates measured at near and far detectors, we 
can extract the oscillation probability between the two detectors that is purely due to the conventional three-neutrino 
oscillation. In addition, fitting the data with a parametrizable flux depletion/excess, such as gaussian-like shape as 
implemented in [25], can also correct the best fit-value. 

In Fig. 5 we show the prior (input) and posterior (output) distributions of the solar mixing angle at JUNO, 
simulated with NuPro [9]. In total, the whole range from 1.8 MeV to 8 MeV is divided into 400 bins equally. Since the 
event rate in each bin is more that 20 events, the statistical fluctuation can be approximated as gaussian distribution. 
The y; 2 -function is defined as 


= V (jV * zJ Ni)2 

“ Ni 


Am? — Am? 


-1 2 r 


J(Am?) 


Am? — Am? 


-1 2 


<5(Am?) 


sin 2 2 9 r — sin 2 2 9 r 


d(sin 2 9 r ) 


sin 2 6 b — sin 2 0 , 


d(sin“ 9 S ) 


1 2 


+ 


(/- l ) 2 

0.01 2 


(3.3) 
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Fig. 4: The event rate of reactor neutrinos to be observed at JUNO. Two values of the solar neutrino mixing angle 
were used, sin 2 0 S = 0.323 (thick lines) and sin 2 0 S = 0.339 (thin lines). 
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Fig. 5: Prior (input) and posterior (output) distributions of the solar mixing angle 0 S at JUNO. 


with input values, 


Am 2 = (7.6 ± 0.2) x 10 -5 eV 2 , Am 2 = (2.4 ± 0.1) x 1CT 3 eV 2 , (3.4a) 

sin 2 2 9 r = 0.089 ± 0.005 , sin 2 9 S = 0.323 ± 0.016 , (3.4b) 

according to the latest global fit [2] and Daya Bay measurement on 9 r [30]. In addition, the flux has an overall 
1% uncertainty parametrized as an overall normalization /. Instead of minimization, the y 2 -function (3.3) is used to 
sample the oscillation parameters, two mass-squared differences Am 2 and Am 2 as well as the two mixing angles 9 r 
and 9 Sl with the Bayesian Nested Sampling algorithm [31]. Since there is no measurement of the two Majorana CP 
phases, they are randomly sampled between 0 and 27 t. The sampled points are then analyzed to get the best fit values, 
uncertainties, and ranges of various functions of oscillation parameters simultaneously. In this way, the maginalization 
is carried out automatically. 

It turns out that JUNO can significantly improve the uncertainty on 9 S from A[sin 2 0 S ] = 0.016 down to 2.4 x 10~ 3 
which is consistent with the rough estimation based on total event numbers. Correspondingly, the uncertainty A (9 S ) 
is as small as 0.15°, roughly 4 x 10~ 3 of 9 S itself. At this precision, we can claim that the solar mixing angle 9 S as well 
as cos 2 9 S are precisely known. In addition, the atmospheric mass-squared difference Am 2 is measured with precision 
A[Am 2 ]/Am 2 ~ 6 x 10~ 3 . In other words, there is almost no contribution from the neutrino mixing sector to the 
global lower limit of as defined in (2.4). This will be quantified in the next Section. 
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Fig. 6 : The 3 a range of (m„) m i n as a function of the lightest neutrino mass, mi for NO and m 3 for 10. All oscillation 
parameters, including the solar and reactor mixing angles, 9 S and 6 ri as well as the atmospheric and solar mass-squared 
differences, Am 2 and Am,, are varied in their corresponding 3er ranges. 


4. Implications for Double Beta Decay 

With the solar mixing angle 9 S and the atmospheric mass-squared difference Am 2 precisely measured by huge statistics 
and excellent energy resolution, respectively, the global lower limit on the effective electron neutrino mass (2.4) is 
actually fixed. In Fig. 6 , we show the predicted as a function of the smallest neutrino mass eigenvalue, in this 

case m 3 , with all other parameters, including the solar and reactor mixing angles, 9 S and 9 r , as well as the atmospheric 
and solar mass-squared differences, AmJ and Am 2 , varied in their corresponding present and future 3<r ranges 2 . At 
any value of m 3 , the uncertainty of the lower limit (m y )^ n is compressed to negligible size. Without JUNO, the value 
of (m l/ )J^ n is unknown within a factor of roughly 1.6 at 3cr confidence level. After JUNO is included, it drops to the 
level of 2 x 10” 3 , which is a huge improvement. It is safe to claim that (m,)^ can be predicted almost precisely. 
For NO, the improvement is also significant, as depicted in Fig. 6 . We also show in Fig. 1 the effective mass versus 
the smallest mass eigenvalue, as well as versus the neutrino mass parameters that are accessible in direct searches and 
cosmology, for both mass orderings. The improvement from the current 3c ranges is obvious. 

Fig. 7 shows the necessary half-lifes to rule out the inverted ordering for the isotopes 48 Ca, 76 Ge, 82 Se, 100 Mo, 
116 Cd, 130 Te, 136 Xe and 150 Nd. We have used a compilation of nuclear matrix elements which is summarized in Table 
1. Also in this plot the improvement from the current range to the future situation is obvious. The shorter error bar 
corresponds to future uncertainty, which solely comes from the nuclear matrix elements. For instance, the range to go 
below (m)J^ n for 76 Ge is currently T®J 2 = (3.7 ... 52.3) x 10 27 yrs, while after JUNO it would be Ty 2 = (5-3 ... 31) x 10 27 
yrs, the range is smaller by a factor 2.4. For 136 Xe the range is currently TW 2 = (1.1... 14.3) x 10 2 ‘ yrs, which would 
improve to T®Y 2 = (1.5... 8.4) x 10 27 yrs. See also Fig. 3 for the necessary half-lifes to rule out (and reach, where 
the experimental improvement is of minor importance) the inverted ordering. For convenience, the ranges of half-lifes 
have been summarized in Table 2. 


5. Conclusion 

The next generation of neutrino oscillation experiments will be able to provide mixing parameter determinations with 
remarkable precision. This has impact on a variety of aspects. Here we have studied inasmuch future reactor neutrino 
experiments can help in planing and interpreting future searches for neutrinoless double beta decay. The minimal 
value of the effective mass in the inverted ordering is currently uncertain by a factor of almost 2 , comparable to the 
nuclear matrix element uncertainty. Facilities like JUNO can determine the minimal value of the effective mass with 
essentially no uncertainty, which fixes the half-life values corresponding to those extreme values, up to the nuclear 
matrix element uncertainties. The total uncertainty is therefore reduced by a factor 2, leaving further improvement 
to the nuclear physics community. 

2 For completeness, we also give in Fig. 6 the corresponding improvement for the minimal value in the normal ordering. 
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Fig. 7: The half-life TjV of neutrinoless double beta decay (0^2 0) to rule out the inverted ordering and its uncertainties. 
If one error bar is displayed, e.g. for SkM, the range is according to the current 3cr uncertainty in the oscillation 
parameters, which will be reduced to a single point after JUNO. If two error bars are displayed, e.g. for IBM-2, the 
nuclear matrix elements have a range, see Table 1. The larger error bar is then the combined uncertainty, from the 
matrix elements and the oscillation parameters in (m)[9 n . 
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QRPA/T [3 1 

QRPA/J [: 1 pnQRPA [39 l 

ISM 1 " 1 ' 1 

IBM-2 [:! - ] 

PHFE 

[35] 

EDF [ 1 

qOv [40] 

20 Ca 

- 

0.541.. 

.0.594 

- 


- 

0.61 

1.75 

X 

(1 

± 

0.16) 


- 


2.37 

23.29 

ItGe 

5.09 

5.157.. 

.5.571 

4.029. . 

.5.355 

5.26 

2.30 

4.68 

X 

(1 

± 

0.16) 


- 


4.60 

2.218 

11 Se 

- 

4.642.. 

.5.018 

2.771.. 

.3.722 

3.73 

2.18 

3.73 

X 

(1 

± 

0.16) 


- 


4.22 

9.537 

'SSMo 

- 

5.402.. 

.5.850 

2.737.. 

.3.931 

3.90 

- 

4.22 

X 

(1 

± 

0.16) 

6.26 

± 

0.63 

5.08 

14.94 

“ICd 

- 

4.040.. 

.4.367 

3.034.. 

.3.935 

4.26 

- 

3.10 

X 

(1 

± 

0.16) 
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4.72 

15.68 

1 i§Te 

1.37 

3.888.. 

.4.373 

2.993.. 

.4.221 

4.00 

2.12 

3.70 

X 

(1 

± 

0.16) 

4.05 

± 

0.49 

5.13 

13.35 

'IfXe 

1.89 

2.177.. 

.2.460 

2.053.. 

.2.802 

2.91 

1.77 
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(1 

± 

0.16) 


- 


4.20 

13.69 

igNd 

2.71 



- 


- 

- 

2.67 
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(1 
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0.16) 

2.83 
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0.42 

1.71 

59.16 


Table 1: Theoretical predictions for the nuclear matrix elements and the phase space factor in units of 10 26 yr 1 eV 2 
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T 1/2 [yrs] 

w/o JUNO 

= (0.0 1 27...0.0 1 98) eV 

with JUNO 

MIL = 0-0166 eV 

NME 

48 Ca 

(1.95...4.69) x 10 27 
(3.74...9.01) x 10 28 

2.76 x 10 27 

5.30 x 10 28 

<a X = 2.37 
<n = 0.55 

76 Ge 

(3.70...8.92) x 10 27 
(2.17...5.23) x 10 28 

5.25 x 10 27 

3.08 x 10 28 

M 3 - ax = 5.57 
<in = 2-30 

82 Se 

(1.06...2.56) x 10 27 
(5.62... 13.5) x 10 27 

1.51 x 10 27 

7.98 x 10 27 

AC ax = 5.02 

M^ n = 2-18 

100 Mo 

(3.59...8.65) x 10 26 
(2.28...5.48) x 10 27 

5.10 x 10 26 

3.23 x 10 27 

AC ax = 6.89 
<in = 2-74 

116 Cel 

(7.29... 17.6) x 10 26 
(2.28...5.48) x 10 27 

1.04 x 10 27 

3.23 x 10 27 

M°" ax = 4.72 
<n = 2-67 

130 Te 

(7.25... 17.5) x 10 26 
(1.02...2.45) x 10 28 

1.03 x 10 27 

1.44 x 10 28 

AC ax = 5.13 

M min = 1-37 

136 Xe 

(1.06...2.54) x 10 27 
(5.94... 14.3) x 10 27 

1.50 x 10 27 

8.43 x 10 27 

M ma X = 4-20 

Cn = 1-77 

150 Nd 

(4.08...9.82) x 10 26 
(1.47...3.55) x 10 27 

5.78 x 10 26 

2.09 x 10 27 

M ma X = 3-25 

M mi n = 1-71 


Table 2: Half-life ranges to rule out the inverted ordering before (after) JUNO. For each isotope, the first row 
corresponds to the range (or value) for the largest matrix element M^ ax , while the second row is the range (value) for 
the smallest matrix element M^ in . Before JUNO (w/o), varies due to the particle physics uncertainty which 

essentially vanishes after JUNO. 
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